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Abstract 

FOCUS results from Dalitz-plot analyses of Df and to tt'^tt^tt'*" are presented. 
The K-matrix formalism is applied to charm decays for the first time, which allows 
us to fully exploit the already existing knowledge coming from light-meson spec- 
troscopy experiments. In particular all the measured dynamics of the >S-wave tttt 
scattering, characterized by broad/overlapping resonances and large non-resonant 
background, can be properly included. This paper studies the extent to which the 
K-matrix approach is able to reproduce the observed Dalitz plot and thus help us 
to understand the underlying dynamics. The results are discussed along with their 
possible implications for the controversial a meson. 

Key words: Amplitude analysis, charm decay, light scalars 
PACS: 



1 Introduction 

Charm-meson decay dynamics has been extensively studied in the last decade. 
The analysis of the three-body final state by fitting Dalitz plots has proved 
to be a powerful tool for investigating effects of resonant substructure, inter- 
ference patterns, and final state interactions in the charm sector. The isobar 
formalism, which has traditionally been applied to charm amplitude analyses. 
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represents the decay amplitude as a sum of relativistic Breit-Wigner propa- 
gators multiplied by form factors plus a term describing the angular distribu- 
tion of the two body decay of each intermediate state of a given spin. Many 
amplitude analyses require detailed knowledge of the light-meson sector. In 
particular, the need to model intermediate scalar particles contributing to the 
charm meson in the decays reported here has caused us to question the validity 
of the Breit-Wigner approximation for the description of the relevant scalar 
resonances [1,2]. Resonances are associated with poles of the S-matrix in the 
complex energy plane. The position of the pole in the complex energy plane 
provides the fundamental, model-independent, process- independent resonance 
description. A simple Breit-Wigner amplitude corresponds to the most ele- 
mentary type of extrapolation from the physical region to an unphysical-sheet 
pole. In the case of a narrow, isolated resonance, there is a close connection 
between the position of the pole on the unphysical sheet and the peak we ob- 
serve in experiments at real values of the energy. However, when a resonance 
is broad and overlaps with other resonances, then this connection is lost. The 
Breit-Wigner parameters measured on the real axis (mass and width) can 
be connected to the pole-positions in the complex energy plane only through 
models of analytic continuation. 

A formalism for studying overlapping and many channel resonances has been 
proposed long ago and is based on the K-matrix [3,4] parametrization. This 
formalism, originating in the context of two-body scattering, can be gener- 
alized to cover the case of production of resonances in more complex reac- 
tions [5], with the assumption that the two-body system in the final state is 
an isolated one and that the two particles do not simultaneously interact with 
the rest of the final state in the production process [4] . The K-matrix approach 
allows us to include the positions of the poles in the complex plane directly 
in our analysis, thus directly incorporating the results from spectroscopy ex- 
periments [6, 7] . In addition, the K-matrix formalism provides a direct way 
of imposing the two-body unitarity constraint which is not explicitly guaran- 
teed in the simple isobar model. Minor unitarity violations are expected for 
narrow, isolated resonances but more severe ones exist for broad, overlapping 
states. The vahdity of the assumed quasi two-body nature of the process of the 
K-matrix approach can only be verified by a direct comparison of the model 
predictions with data. In particular, the failure to reproduce three-body-decay 
features would be a strong indication of the presence of the neglected three- 
body effects. 



2 Candidate selection 

The FOCUS detector is a large aperture, fixed-target spectrometer with ex- 
cellent vertexing and particle identification capabilities. We have chosen cuts 
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designed to minimize non-charm background as well as reflection backgrounds 
from misidentified charm decays. The three-pion final states are obtained us- 
ing a candidate driven vertex algorithm. A decay vertex is formed from three 
reconstructed charged tracks. The momentum of the D candidate is used to 
intersect other reconstructed tracks to form a production vertex. The confi- 
dence levels (C.L.) of each vertex is required to exceed 1%. After the vertex 
finder algorithm, the variable i, which is the separation of the primary and 
secondary vertex, and its associated error cr^ are calculated. We reduce back- 
grounds by requiring i/ai > 10 and 7 for the D"*" and D+, respectively. The 
two vertices are also required to satisfy isolation conditions. The primary ver- 
tex isolation cut requires that a track assigned to the decay vertex has a C.L. 
less than 1 % to be included in the primary vertex. The secondary vertex iso- 
lation cut requires that all remaining tracks not assigned to the primary and 
secondary vertex have a C.L. smaller than 0.1 % to form a vertex with the D 
candidate daughters. The decay vertex is required to be 3 cr outside of the tar- 
get material to reduce the background due to hadronic re-interactions in the 
material. A cut on the negative log likelihood of the Cerenkov hypothesis [8] 
of Wk — W^TT > is required for each pion. A tighter cut of Wk — W^^ > 5 is re- 
quired on opposite-sign pion in the decay in order to remove the K~t:^t:^ 
reflection contribution to the low-mass sideband. We further require that 
all three pions satisfy a loose pion-consistency cut of AVF = W.„ — W^an < 3 
where IVmin is the negative log likelihood of the most favored Cerenkov hy- 
pothesis. The vertex isolation requirement nearly eliminates D*^ — > D^tt~^ 
contamination. The samples selected according to these requirements (Fig. 1) 
consist of 1475 ± 50 and 1527 ± 51 signal events for the Df and respec- 
tively. The Dalitz plot analyses are performed on events within a the nominal 
or mass (Fig. 2). 

3 The decay amplitude 

The decay amplitude of the D meson into the three-pion flnal state is written 
as: 

A{D) = aoe''° + F, + ^ a,e^^'S(a6c|r,), (1) 

i 

where the first term represents the direct non-resonant three-body amplitude 
contribution, Fi is the contribution of 5'-wave states and the sum is over the 
contributions from the intermediate two-body non-scalar resonances. B{abc\ri) 
are the usual Breit-Wigner terms of the traditional isobar model, whose ex- 
plicit forms are given in [9]. Fi is written in the context of the K-matrix 
approach which we will discuss shortly. The coefficients and phases of the 
A{D) amplitude are all relative to a free parameter of the Fi amplitude, 
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Figure 1. Signal and side-band regions of the three-pion invariant-mass distribution 
for a) and b) Dalitz-plot analysis respectively. 
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Figure 2. a) Df and b) Dalitz plots. 

whose modulus and phase are fixed to 1 and respectively (see below). The 
squared modulus of this amplitude gives the probability density in the three 
pion Dalitz plot. 

In general, the decay of a D meson into three pions via a resonance involves 
the production of a {IJ)^^ state with an accompanying pion. It is believed 
that, while the usual Breit-Wigner (BW) approximation is suitable for states 
with J > 0, since they are characterized by relatively narrow and isolated 
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resonances, the treatment of 5'-wave states requires a more general formalism 
to account for non-trivial dynamics due to the presence of broad and overlap- 
ping resonances [1,2]. For J = 0, only states with even isospin and positive P 
and C are allowed to strongly couple to tt'^'tt". We limit ourselves to isoscalar 
S-w&ye states, (00)'^"'", since 1 — 2 must involve at least two qq pairs and no 
four-quark states with 1 = 2 are known. At the mass scales relevant to this 
analysis, the decay of a charm particle into a (00)"^+ state with an accom- 
panying pion consists of five channels / = 1 ... 5 where 1 = tttt, 2 = KK, 
3 = multi-meson states (four-pion state mainly at y/s < 1.6 GeV) i = rjr] and 
5 = 7777'. The amplitude for the particular channel (00)^^7r can be written in 
the context of the K-matrix formalism as 

Fi = {I- iKp)i^P, , (2) 



where / is the identity matrix, K is the K-matrix describing the isoscalar 
(S-wave scattering process, p is the phase-space matrix for the five channels, 
and P is the "initial" production vector into the five channels. In this picture, 
the production process can be viewed as consisting of an initial preparation 
of several states, which are then propagated by the (/ — iKp)~^ term into the 
final one. Only the Fi amplitude is present in the isosinglet 5'-wave term since 
we are describing the dipion channel. 

We require a reliable K-matrix parametrization of (00)++-wave scattering. 
To our knowledge the only self-consistent description of the ^-wave isoscalar 
scattering is given by the K-matrix representation of Anisovich and Sarantsev 
in reference [7] obtained through a global fit of the available scattering data 
from TTTT threshold up to 1900 MeV. Their K-matrix parametrization is: 



The 51}" is the coupling constant of the K-matrix pole rria to the i meson 
channel; the parameters ffj^^^ and s^^^^ describe a slowly varying part (which 
we will call SVP) of the K-matrix elements; the factor i^-^Aml/2){i-SAo) j 
to suppress false kinematical singularity in the physical region near the tttt 
threshold (Adler zero). The parameter values used in this paper are listed 
in Table 1, which was provided by the authors of reference [7]. Note that 
K-matrix representation is by definition real and symmetric. The K-matrix 
values of Table 1 generate a physical T-matrix, T = {I — iK ■ p)~^K, which 
describes the scattering in the (00)++- wave with five poles, whose masses, 
half- widths, and couphngs are listed in Table 2. 

The /o series reported in Table 2 differs somewhat from that reported by 
the PDG [10] group. In addition to the /o(980) and /o(1500) poles which 
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Table 1 

K-matrix parameters. Masses and coupling constants are in GeV. Only the i = 1 
fij terms are reported here since they are the only values relevant to the three-pion 
decay. 

"f^a S'tttt 9kK 947r S'7777 9rjr)' 

0.65100 0.24844 -0.52523 0.00000 -0.38878 -0.36397 

1.20720 0.91779 0.55427 0.00000 0.38705 0.29448 

1.56122 0.37024 0.23591 0.62605 0.18409 0.18923 

1.21257 0.34501 0.39642 0.97644 0.19746 0.00357 

1.81746 0.15770 -0.17915 -0.90100 -0.00931 0.20689 

„scatt j^scatt j-scatt fscatt /-scatt /-scatt 

*0 /ll /l2 /l3 /l4 Jib 

-3.30564 0.26681 0.16583 -0.19840 0.32808 0.31193 

SA SAO 

1.0 -0.2 



Table 2 



T-matrix poles and tt+tt" scattering coupling-constants; phases are in degrees. The 
/o(980) pole position is found by the authors of [7] on the second Riemann sheet, 
thus the closest physical region is the beginning of KK threshold. 



T-matrix pole 
( m,r/2 ) GeV 


9kkI 91^1^ 




9r)ri 1 91^1^ 


9rp]' 1 9tvtv 


/o(980) 
(1.019,0.038) 


1.3970 e*^^-^ 


0.3572 e*^^-^ 


1.1660 e*^^-5 


0.9662 e*^^° 


/o(1300) 
(1.306,0.170) 


0.2579 e-*i6-6 


2.1960 e-'i'^s'^ 


0.3504 e^23-2 


0.5547 e*i6-2 


/o(1200 - 1600) 
(1.470,0.960) 


1.1140 e-^°-^ 


2.2200 e-^6-^ 


0.569 e^^'^'^ 


0.2309 e^5^-6 


/o(1500) 
(1.488,0.058) 


0.5460 e-'^-^ 


1.9790 e^^^-^ 


0.4083 e^^^'' 


0.4692 e^^^-6 


/o(1750) 
(1.746, 0.160) 


0.1338 e*3^-^ 


1.3680 e'134.8 


0.2979 e'25 i 


0.5843 e"'°-^ 



also appear in the PDG classification, three other poles are present, /o(1300), 
/o(1200 — 1600) and /o(1750), in contrast with only two poles listed by the 
PDG, /o(1370) and /o(1710). The five-pole series used here is able to con- 
sistently reproduce the available ^-wave isoscalar data in the energy range 
relevant for this analysis. The decay amphtude for the D meson into the 
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three-pion final state, where tt+tt are in a {IJ^^ — 00+''")-wave is then 



w -(T n R-n^-i J ^"^i I /prod 1 GeV - si 

a '"'a * S — Sq 



^ {s- SAml/2){l- Sao) 
(s - Sao) 

where (3a is the couphng to the nia pole in the 'initial' production process, 
and Sq™'^ are the P-vector SVP parameters. (3^ and ffj°'^ are in general complex 
numbers [4] . The phase space matrix elements for the two pseudoscalar-particle 
states are: 



\ 



mu + m2i)^\ {mu-m2 




TTTT, KK, rjrj, rjr]') . (5) 



The normalization is such that — 1 as s — cxd. The expression for the 
multi-meson state phase space can be found in reference [7]. 

We note that the P-vector poles have to be the same as those of the K-matrix 
in order to cancel out infinities in the final amplitude as each pole is realized. 
The P-vector SVP parametrization is chosen in complete analogy with that 
used for the K-matrix. The need for the Adler-zero term, not a-priori required 
in the P-vector, will be investigated by studying its effect on the quality of the 
fit to our data. The K-matrix parameters are fixed to the values of Table 1 
in our Dalitz plot fits. The free parameters are the P-vector parameters {(3a, 
fij°'^ and Sq™'^), and the coefficients and phases of Eq. 1 (oq, and ^o, ^i)- All 
amplitudes are referenced to (3i which is fixed at 1. The P-vector Adler-zero 
parameters have been chosen to be identical to those of the K-matrix, sa and 
Sao- Because of the limited allowed range for these parameters, we do not 
expect our results to critically depend on this particular choice. 



4 The likelihood function and fitting procedure 



The probability density function is corrected for geometrical acceptance and 
reconstruction efficiency. We find that finite-mass resolution effects are negligi- 
ble. The shape of the background in the signal region is parametrized through 
a polynomial fit to the Dalitz plot of mass sidebands ^ . The number of back- 

^ For this analysis the sideband between the two signal peaks begins at 3 a from 
the D'^ peak and ends at 3 cr from the Df peak where the cr's are the r.m.s. widths 
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ground events expected in the signal region is estimated through fits to the 
7r"'"7r~7r"'" mass spectrum. All background parameters are included as additional 
fit parameters and tied to the results of the sideband fits through the inclusion 
in the likelihood of a penalty term derived from the covariance matrix of the 
sideband fit. The contamination in the left sideband from — > K~'k'^'k'^ ^ 
where is misidcntificd as tt", is reduced to a negligible level (3.5 % of the 
total events in the sideband) using the tight Ccrcnkov cut. Background from 
the decay, — >• r^'TT"*" with rj' ^ fq and p — >• 71^71^ , is expected in the 
signal and sideband regions. It is included by adding a p-BW component 
in the background parametrization. The and samples are fitted with 
likelihood functions C consisting of signal and background probability den- 
sities. Checks for fitting procedure are made using Monte Carlo techniques 
and all biases are found to be small compared to the statistical errors. The 
systematic errors on our results are evaluated by comparing their values in 
disjoint samples corresponding to different experimental running conditions, 
different kincmatical regions, such as low versus high D momenta, and particle 
versus anti-particle. A split sample systematic error was added in quadrature 
to the existing statistical error to make the split sample estimates consistent 
to within 1 cr if necessary. The assumption that the shape of the background 
in the sideband is a good representation of the background in the signal re- 
gion could potentially constitute another source of systematic error. We study 
this effect by varying the polynomial function degree and adding/removing 
the Breit-Wigner terms, which are introduced to take into account any feed- 
through from resonances in the background, and computing the r.m.s. of the 
different results, which is added in quadrature to form the total experimental 
systematic error. 



5 Results for the £>+ — > tt+tt 7r+ decay 

We recall that the physical parameters of our fit are P-vector parameters: (3^, 
ff'j°'^j along with Sq™'^, and the coefficients and phases of Eq. 1, aQ,ai and 
do,di. The K-matrix parameters are fixed to the values given in Table 1. The 
general procedure, adopted for all the fits reported here, consists of several 
successive steps in order to eliminate contributions whose effects on our fit are 
marginal. We initially consider all the well established, non-scalar resonances 
decaying to tt+tt" with a sizeable branching ratio. Contributions are removed 
if their amplitude coefficients, of Eq. 1, arc less than 2 a significant and 
the fit confidence level increases due to the decreased number of degrees of 
freedom in the fit. The P-vector initial form includes the complete set of K- 

of the two measured mass peaks. The left sideband for covers the —5 a to — 3(T 
region from the peak, while the right sideband for Df, the 3 a" through 6 a 
region from the Df mass peak. 
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matrix poles and slowly varying function (SVP) as given in reference [7]; (3a 
as well as the ffj°'^ terms of Eq. 4 arc removed with the same criteria. The 
fit confidence levels (C.L.) are evaluated with a estimator over a Dalitz 
plot with bin size adaptively chosen to maintain a minimum number of events 
in each bin. Once the minimal set of parameters is reached, addition of each 
single contribution previously eliminated is reinstated to verify that the C.L. 
does not improve. 

Table 3 shows the P-vector composition from our final fit results on the Df — > 
7r+7r~7r+ Dahtz plot. The fifth K-matrix pole and the second SVP contribution 

Table 3 



P-vector parameters of the Df fit. 


P-vector parameters 


modulus 


phase (deg) 


Pi 


1 (fixed) 


(fixed) 




0.974 ±0.102 


207.0 ± 5.9 


03 


0.755 ± 0.196 


85.7 ± 19.3 


04 


1.760 ±0.371 


119.1 ± 19.9 


/•prod 

Ju 


0.863 ±0.153 


202.3 ± 9.4 


j^prod 
^13 


3.768 ± 0.746 


140.3 ± 16.7 


(•prod 
JU 


2.179 ± 0.663 


106.6 ± 24.4 


/•prod 
/l5 


2.452 ± 0.684 


236.5 ± 13.7 



were eliminated. The inclusion of an Adler zero term did not improve our fit 
quality and was removed. The quoted results were obtained with s'q°'^ — —3 
GeV^, but they were insensitive to any choice in the range —5 GeV^ < Sq^°'^ < 
— 1 GeV^ - typical parameter values for the SVP. 

The resulting fit fractions ^ , phases and amplitude coefficients are quoted 
in Table 4. We note that both the three-body non- resonant and p°(770)7r"'" 
components were not required by the fit. We represent the entire S-wave con- 
tribution by a single fit fraction since, as previously discussed, one cannot dis- 
tinguish the different resonance or SVP S'-wave contributions on the real axis. 
The couplings to T-matrix physical poles, reported in Table 5, are computed 
by continuing the amplitude Fi{s) into the complex s- plane to the position 
of the poles and evaluating the pole residues ^ . The Dahtz projections 

^ The quoted fit fractions are defined as the ratio between the intensity for a single 
amplitude integrated over the Dalitz plot and that of the total amplitude with all 

the modes and interferences present. 

^ The coupling to the /o(1200— 1600) pole has to be looked at with a certain caution 
because of the intrinsic limitations of the approximation used for the 4-body phase- 
space when extrapolated very deeply into the complex plane. 
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of our data are shown in Fig. 3 superimposed with our final fit projections. 
Figure 4 shows the corresponding adaptive binning scheme used to obtain the 
fit confidence level. 

Table 4 



Fit results from the K-matrix model for Df. 


decay channel 


fit fraction (%) 


phase (deg) 


amplitude coefficient 


(5- wave) 7r+ 


87.04 ±5.60 ±4.17 


(fixed) 


1 (fixed) 


/2(1270) 7r+ 


9.74 ±4.49 ±2.63 


168.0 ±18.7 ±2.5 


0.165 ±0.033 ±0.032 


/(1450)7r+ 


6.56 ±3.43 ±3.31 


234.9 ± 19.5 ± 13.3 


0.136 ±0.030 ±0.035 


Fit C.L 


3.0% 







Table 5 

7r+7r~7r"'" production coupling constants for the five T-matrix poles, which 
comprise the S-wave component in Table 4, referenced to the /o(980). Phases are 
in degrees. 



T-matrix pole 


(m,r/2) (GeV) 


D'^ (relative) coupling constant 


/o(980) 


(1.019,0.038) 


1 e*° (fixed) 


/o(1300) 


(1.306,0.170) 


(0.43 ± 0.04) e^(-i63.8±4.9) 


/o(1200 - 1600) 


(1.470,0.960) 


(4.90 ± 0.08) e*( 80-9±i-06) 


/o(1500) 


(1.488,0.058) 


(0.51 ± 0.02) e^( 83.i±3.03) 


/o(1750) 


(1.746,0.160) 


(0.82 ± 0.02) e*(-127.9±2.25) 




Figure 3. Df Dalitz-plot projections with our final fit superimposed. The back- 
ground shape under the signal is also shown. 
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Figure 4. Df adaptive binning Dalitz-plot for the K-matrix model. 
6 Results for the £)+ — > 7r+7r~7r+ decay 

The — > 7r"'"7r~7r"'" Dahtz plot shows an excess of events at low ti^ti~ mass, 
which cannot be explained in the context of the simple isobar model with 
the usual mixture of well established resonances along with a constant, non- 
resonant amplitude. A new scalar resonance, the a(600), has been previously 
proposed [11] to describe this excess. However we know that complex structure 
can be generated by the interplay among the S-wave resonances and the un- 
derlying non-resonant S-wave component that cannot be properly described 
in the context of a simple isobar model. It is therefore interesting to study 
this channel with the present formalism, which embeds all our experimental 
knowledge about the S'-wave vr+vr" scattering dynamics. 

With the same procedure based on statistical significance and fit confidence 
level used in the analysis, we obtained the final set of P-vector parameters 
that is reported in Table 6. The last two poles and the last three SVP terms 

Table 6 



P-vector parameters 


of the D+ fit. 




P-vector parameters modulus 


phase (deg) 


Pi 


1 (fixed) 


(fixed) 


P2 


2.471 ±0.431 


82.5 ± 10.3 


P3 


1.084 ± 0.386 


102.8 ±23.5 


(■prod 

Jn 


2.565 ±0.737 


155.4 ± 18.3 


(.prod 
/l2 


6.312 ±0.967 


-160.0 ±8.7 
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were eliminated. The Sq™'^ value is measured to be s'q°'^ — {—l.Otl'i) GeV^ . 
The fit did not require an Adler-zero term. 

Beside the S-wavc component, the decay appears to be dominated by the 
p°(770) plus a /2(1270) component. The p*^(1450) was always found to have 
less than 2 a significance and was therefore dropped from the final fit. In 
analogy with the Df, the direct three-body non- resonant component was not 
necessary since the SVP of the S'-wave could reproduce the entire non-resonant 
portion of the Dalitz plot. The complete fit results are reported in Table 7. 
The resulting production couphng constants are reported in Table 8. The 

Table 7 



Fit results from the K-matrix model fit for D~^. 



decay channel 


fit fraction (%) 


phase (deg) amplitude coefficient 


(S'-wave) 7r+ 


56.00 ± 3.24 ± 2.08 


(fixed) 1 (fixed) 


/2(1270)7r+ 


11.74 ± 1.90 ±0.23 


-47.5 ± 18.7 ± 11.7 1.147 ± 0.291 ± 0.047 


p0(770)7r+ 


30.82 ±3.14 ±2.29 


-139.4 ± 16.5 ± 9.9 1.858 ± 0.505 ± 0.033 


Fit C.L. 


7.7% 




Table 8 

D'^ — 7r'^7r~7r"'" production coupling constants for the five T-matrix poles, which 
comprise the S'-wave component in Table 7, referenced to the /o(980). Phases are 
in degrees. 


T-matrix pole 


(m,r/2) (GeV) 


(relative) coupling constant 


/o(980) 


(1.019,0.038) 


1 e^° (fixed) 


/o(1300) 


(1.306,0.170) 


(0.67 ±0.03) e^(-67.9±3.o) 


/o(1200 - 1600) 


(1.470,0.960) 


(1.70 ±0.17) e'(-i25.5±i.7) 


/o(1500) 


(1.488,0.058) 


(0.63 ± 0.02) e'(-i42.2±2.2) 


/o(1750) 


(1.746,0.160) 


(0.36 ± 0.02) e^(-i35-0±2.9) 



Dalitz projections are shown in Fig. 5 and the corresponding adaptive binning 
scheme is shown in Fig. 6. The most interesting feature of these results is 
the fact that the better treatment of the S'-wave contribution provided by the 
K-matrix model can reproduce the low-mass tt+tt" structure of the D'^ Dalitz 
plot. This suggests that any o"-like object in the D decay should be consistent 
with the same cr-like object measured in the tt+tt^ scattering. We believe 
that additional studies with higher statistics will be required to completely 
understand the a puzzle. 
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Low ^+ ^- mass projection (GeV^ Ic^) High ^+ ^- mass projection (GeV^ Ic^) 



Figure 5. Dalitz-plot projections with our final fit superimposed. The back- 
ground shape under the signal is also shown. 



= 48 .673 

Free Pai= 17 

DOF = 36 

: 7.72E-02 



X^/dof = 1.352 



Bins = 53 



Z <1 
1>X^<2 




Figure 6. adaptive binning Dalitz-plot for the K-matrix fit. 



7 £>+ and Z>+ final results 

s 



The K-matrix parameters used in this analysis correspond to the best solution 
provided by the authors of reference [7]. Several solutions with slightly different 
parametrizations for the 4tt phase-space and for the K-matrix background 
terms were presented in the same paper. We evaluate the systematic error 
due to solution choice by computing the r.m.s. of the fit fractions and phases 
obtained using the different solutions. The final results, including this last 
systematic error, are presented in Table 9. 
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Table 9 

Final results on Df and tt~^tt~tt^ fit fractions and phases. Beside the first 
reported error, which is statistical, two systematic errors are quoted. The first one 
is from the measurement systematics and the second one is due to the particular 
solution chosen for the K-matrix poles and backgrounds. 



decay channel 


fit fraction (%) 


phase (deg) 


(S-wave) tt'^ 


87.04 ±5.60 ±4.17 ±1.34 


(fixed) 


/2(1270)7r+ 


9.74 ± 4.49 ± 2.63 ± 1.32 


168.0 ±18.7 ±2.5 ±21.7 


/90(1450)7r+ 


6.56 ±3.43 ±3.31 ±2.90 


234.9 ± 19.5 ± 13.3 ± 24.9 


D+ 


decay channel 


fit fraction (%) 


phase (deg) 


(S-wave) TT"*" 


56.00 ± 3.24 ± 2.08 ± 0.50 


(fixed) 


/2(1270)7r+ 


11.74 ± 1.90 ±0.23 ±0.18 


-47.5 ±18.7 ±11.7 ±5.3 




30. 82 ±3.11 ±2. 29 ±0.17 


- 139. 1± If). 5 ±9. 9 ±5.0 



8 Conclusions 



The K-matrix formalism has been applied for the first time to the charm 
sector in our Dalitz plot analyses of the D'^ and — > 7r+7r~7r+ final states. 
The results are extremely encouraging since the same K-matrix description 
gives a coherent picture of both two-body scattering measurements in light- 
quark experiments as well as charm meson decay. This result was not obvious 
beforehand. Furthermore, the same model is able to reproduce features of 
the — > 7r"'"7r~7r+ Dalitz plot that otherwise would require an ad hoc a 
resonance. In addition, the non-resonant component of each decay seems to 
be described by known two-body (S-wave dynamics without the need to include 
constant amplitude contributions. 

The K-matrix treatment of the S-wave component of the decay amplitude 
allows for a direct interpretation of the decay mechanism in terms of the five 
virtual channels considered: tttt, KK, rjrj, rjrj' and 47r. By inserting KK~^ in 
the decay amplitude, F, 

F^{I- iKp)-^P = (/ - iKpY^KK-^P = TK-^P = TQ (6) 

we can view the decay as consisting of an initial production of the five virtual 
states which then scatter via the physical T into the final state. The Q-vector 
contains the production amplitude of each virtual channel in the decay. Fig- 
ure 7 shows the ratio of the moduli of the Q-vector amplitudes with respect 
to the TTTT modulus for the Df S-w&ve. The last plot in Fig. 7 represents 
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the normalizing tttt modulus. The two peaks of the ratios correspond to the 
two dips of the tttt normalizing modulus, while the two peaks due to the K- 
matrix singularities, visible in the normalization plot, cancel out in the ratios. 
Figure 8 shows the analogous plots for the D'^ S-wme decay. The resulting 




Figure 7. Measured Df Q-vector components: the first four plots are the ratios of 
moduU of the Q-vector amplitudes with respect to the tttt modulus; the last plot is 
the normalizing tttt modulus. 

picture, for both and D'^ decay, is that the S'-wave decay is dominated by 
an initial production of rjr], rjr]' and KK states. Dipion production is always 
much smaller. This suggests that in both cases the 5'-wave decay amplitude 
primarily arises from a ss contribution such as that produced by the Cabibbo 
favoured weak diagram for the Df and one of the two possible singly Cabibbo 
suppressed diagrams for the D~^. For the D'^, the ss contribution competes 
with a dd contribution. That the /o(980) appears as a peak in the tttt mass 
distribution in decay, as it does in Dg decay, shows that for the 5- wave 
component the ss contribution dominates [2] . Comparing the relative S'-wave 
fit fractions that we observe for Df and D'^ reinforces this picture. The S'- 
wave decay fraction for the D'^ (87%) is larger than that for the (56%). 
Rather than coupling to an S'-wave dipion, the dd piece prefers to couple to a 
vector state like p°(770) that alone accounts for ~ 30 % of decay. 

This interpretation also bears on the role of the annihilation diagram in the 

7r+7r^7r+ decay. We believe that Fig. 7 suggests that the S-wave an- 
nihilation contribution is negligible over much of the dipion mass spectrum. 
It might be interesting to search for annihilation contributions in higher spin 
channels, such as p°(1450)7r and /2(1270)7r. 
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1 2 3 

|KK| / s(GevV) 



1 2 3 

|47t| / iJtJtl s(GeVV) 




7171 Q-vector modulus s(Gev^/c'') 



Figure 8. Measured Q-vector components: the first four plots are the ratios of 
moduH of the Q-vector amplitudes with respect to the tttt modulus; the last plot is 
the normalizing tttt modulus. 
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